Comparative genomics incorporating translocation renal cell carcinoma mouse model reveals molecular mechanisms of tumorigenesis.

Translocation renal cell carcinoma (tRCC) most commonly involves an ASPSCR1-TFE3 fusion, but molecular mechanisms remain elusive and animal models are lacking. Here, we show that human ASPSCR1-TFE3 driven by Pax8-Cre (a credentialed clear cell RCC driver) disrupted nephrogenesis and glomerular development, causing neonatal death, while the clear cell RCC failed driver, Sglt2-Cre, induced aggressive tRCC (as well as alveolar soft part sarcoma) with complete penetrance and short latency. However, in both contexts, ASPSCR1-TFE3 led to characteristic morphological cellular changes, loss of epithelial markers, and an epithelial-mesenchymal transition. Electron microscopy of tRCC tumors showed lysosome expansion, and functional studies revealed simultaneous activation of autophagy and mTORC1 pathways. Comparative genomic analyses encompassing an institutional human tRCC cohort (including a hitherto unreported SFPQ-TFEB fusion) and a variety of tumorgraft models (ASPSCR1-TFE3, PRCC-TFE3, SFPQ-TFE3, RBM10-TFE3, and MALAT1-TFEB) disclosed significant convergence in canonical pathways (cell cycle, lysosome, and mTORC1) and less established pathways such as Myc, E2F, and inflammation (IL-6/JAK/STAT3, interferon-γ, TLR signaling, systemic lupus, etc.). Therapeutic trials (adjusted for human drug exposures) showed antitumor activity of cabozantinib. Overall, this study provides insight into MiT/TFE-driven tumorigenesis, including the cell of origin, and characterizes diverse mouse models available for research.

Population-wide gene disruption in the murine lung epithelium via AAV-mediated delivery of CRISPR-Cas9 components.

With the aim of expediting drug target discovery and validation for respiratory diseases, we developed an optimized method for in situ somatic gene disruption in murine lung epithelial cells via AAV6-mediated CRISPR-Cas9 delivery. Efficient gene editing was observed in lung type II alveolar epithelial cells and distal airway cells following assessment of single- or dual-guide AAV vector formats, Cas9 variants, and a sequential dosing strategy with combinatorial guide RNA expression cassettes. In particular, we were able to demonstrate population-wide gene disruption within distinct epithelial cell types for separate targets in Cas9 transgenic animals, with minimal to no associated inflammation. We also observed and characterized AAV vector integration events that occurred within directed double-stranded DNA break sites in lung cells, highlighting a complicating factor with AAV-mediated delivery of DNA nucleases. Taken together, we demonstrate a uniquely effective approach for somatic engineering of the murine lung, which will greatly facilitate the modeling of disease and therapeutic intervention.

Determinants of renal cell carcinoma invasion and metastatic competence.

Metastasis is the principal cause of cancer related deaths. Tumor invasion is essential for metastatic spread. However, determinants of invasion are poorly understood. We addressed this knowledge gap by leveraging a unique attribute of kidney cancer. Renal tumors invade into large vessels forming tumor thrombi (TT) that migrate extending sometimes into the heart. Over a decade, we prospectively enrolled 83 ethnically-diverse patients undergoing surgical resection for grossly invasive tumors at UT Southwestern Kidney Cancer Program. In this study, we perform comprehensive histological analyses, integrate multi-region genomic studies, generate in vivo models, and execute functional studies to define tumor invasion and metastatic competence. We find that invasion is not always associated with the most aggressive clone. Driven by immediate early genes, invasion appears to be an opportunistic trait attained by subclones with diverse oncogenomic status in geospatial proximity to vasculature. We show that not all invasive tumors metastasize and identify determinants of metastatic competency. TT associated with metastases are characterized by higher grade, mTOR activation and a particular immune contexture. Moreover, TT grade is a better predictor of metastasis than overall tumor grade, which may have implications for clinical practice.

Single-molecule nanopore sequencing reveals extreme target copy number heterogeneity in arylomycin-resistant mutants.

Tandem gene amplification is a frequent and dynamic source of antibiotic resistance in bacteria. Ongoing expansions and contractions of repeat arrays during population growth are expected to manifest as cell-to-cell differences in copy number (CN). As a result, a clonal bacterial culture could comprise subpopulations of cells with different levels of antibiotic sensitivity that result from variable gene dosage. Despite the high potential for misclassification of heterogenous cell populations as either antibiotic-susceptible or fully resistant in clinical settings, and the concomitant risk of inappropriate treatment, CN distribution among cells has defied analysis. Here, we use the MinION single-molecule nanopore sequencer to uncover CN heterogeneity in clonal populations of Escherichia coli and Acinetobacter baumannii grown from single cells isolated while selecting for resistance to an optimized arylomycin, a member of a recently discovered class of Gram-negative antibiotic. We found that gene amplification of the arylomycin target, bacterial type I signal peptidase LepB, is a mechanism of unstable arylomycin resistance and demonstrate in E. coli that amplification instability is independent of RecA. This instability drives the emergence of a nonuniform distribution of lepB CN among cells with a range of 1 to at least 50 copies of lepB identified in a single clonal population. In sum, this remarkable heterogeneity, and the evolutionary plasticity it fuels, illustrates how gene amplification can enable bacterial populations to respond rapidly to novel antibiotics. This study establishes a rationale for further nanopore-sequencing studies of heterogeneous cell populations to uncover CN variability at single-molecule resolution.

Embryonic lethality and defective mammary gland development of activator-function impaired conditional knock-in Erbb3 V943R mice.

ERBB3 is a pseudokinase domain-containing member of the ERBB family of receptor tyrosine kinases (RTKs). Following ligand binding, ERBB receptors homo- or hetero-dimerize, leading to a head-to-tail arrangement of the intracellular kinase domains, where the "receiver" kinase domain of one ERBB is activated by the "activator" domain of the other ERBB in the dimer. In ERBB3, a conserved valine at codon 943 (V943) in the kinase C-terminal domain has been shown to be important for its function as an "activator" kinase in vitro. Here we report a knock-in mouse model where we have modified the endogenous Erbb3 allele to allow for tissue-specific conditional expression of Erbb3 V943R (Erbb3 CKI-V943R ). Additionally, we generated an Erbb3 D850N (Erbb3 CKI-D850N ) conditional knock-in mouse model where the conserved aspartate in the DFG motif of the pseudokinase domain was mutated to abolish any potential residual kinase activity. While Erbb3 D850N/D850N animals developed normally, homozygous Erbb3 V943R/V943R expression during development resulted in embryonic lethality. Further, tissue specific expression of Erbb3 V943R/V943R in the mammary gland epithelium following its activation using MMTV-Cre resulted in delayed elongation of the ductal network during puberty. Single-cell RNA-seq analysis of Erbb3 V943R/V943R mammary glands showed a reduction in a specific subset of fibrinogen-producing luminal epithelial cells.

Accurate assembly of the olive baboon (Papio anubis) genome using long-read and Hi-C data.

BACKGROUND: Baboons are a widely used nonhuman primate model for biomedical, evolutionary, and basic genetics research. Despite this importance, the genomic resources for baboons are limited. In particular, the current baboon reference genome Panu_3.0 is a highly fragmented, reference-guided (i.e., not fully de novo) assembly, and its poor quality inhibits our ability to conduct downstream genomic analyses. FINDINGS: Here we present a de novo genome assembly of the olive baboon (Papio anubis) that uses data from several recently developed single-molecule technologies. Our assembly, Panubis1.0, has an N50 contig size of ∼1.46 Mb (as opposed to 139 kb for Panu_3.0) and has single scaffolds that span each of the 20 autosomes and the X chromosome. CONCLUSIONS: We highlight multiple lines of evidence (including Bionano Genomics data, pedigree linkage information, and linkage disequilibrium data) suggesting that there are several large assembly errors in Panu_3.0, which have been corrected in Panubis1.0.

Elevated Serum Amino Acids Induce a Subpopulation of Alpha Cells to Initiate Pancreatic Neuroendocrine Tumor Formation.

The cellular origin of sporadic pancreatic neuroendocrine tumors (PNETs) is obscure. Hormone expression suggests that these tumors arise from glucagon-producing alpha cells or insulin-producing ß cells, but instability in hormone expression prevents linage determination. We utilize loss of hepatic glucagon receptor (GCGR) signaling to drive alpha cell hyperproliferation and tumor formation to identify a cell of origin and dissect mechanisms that drive progression. Using a combination of genetically engineered Gcgr knockout mice and GCGR-inhibiting antibodies, we show that elevated plasma amino acids drive the appearance of a proliferative population of SLC38A5+ embryonic progenitor-like alpha cells in mice. Further, we characterize tumors from patients with rare bi-allelic germline GCGR loss-of-function variants and find prominent tumor-cell-associated expression of the SLC38A5 paralog SLC7A8 as well as markers of active mTOR signaling. Thus, progenitor cells arise from adult alpha cells in response to metabolic signals and, when inductive signals are chronically present, drive tumor initiation.

Integrated genomic analysis reveals mutated ELF3 as a potential gallbladder cancer vaccine candidate.

Gallbladder cancer (GBC) is an aggressive gastrointestinal malignancy with no approved targeted therapy. Here, we analyze exomes (n = 160), transcriptomes (n = 115), and low pass whole genomes (n = 146) from 167 gallbladder cancers (GBCs) from patients in Korea, India and Chile. In addition, we also sequence samples from 39 GBC high-risk patients and detect evidence of early cancer-related genomic lesions. Among the several significantly mutated genes not previously linked to GBC are ETS domain genes ELF3 and EHF, CTNNB1, APC, NSD1, KAT8, STK11 and NFE2L2. A majority of ELF3 alterations are frame-shift mutations that result in several cancer-specific neoantigens that activate T-cells indicating that they are cancer vaccine candidates. In addition, we identify recurrent alterations in KEAP1/NFE2L2 and WNT pathway in GBC. Taken together, these define multiple targetable therapeutic interventions opportunities for GBC treatment and management.

Atypical presentation and pathogenesis of a macaque lymphocryptoviral-associated B-cell lymphoma in a cynomolgus monkey.

We report the unique pathogenesis and presentation of a rapidly progressive B-cell lymphoma in a 3-year-old female cynomolgus monkey on day 50 of a 13-week toxicity study. Clinical pathology evaluation revealed a marked leukocytosis with bicytopenia. A serum protein electrophoresis was consistent with monoclonal gammopathy. The architecture of the lymph node, spleen, and thymus were variably effaced by neoplastic cells, which also infiltrated other tissues. Immunohistochemistry of the affected tissues confirmed a predominant population of CD20+, CD79a+, CD3-, CD68-, and CD34-neoplastic cells. The full data best support a diagnosis of Stage V lymphoma. Nextgen sequencing and negative prestudy serology results suggested a recent infection by macaque lymphocryptovirus (mLCV) with a unique transcriptional profile comparable with a rarely observed direct LCV infection model. This infection model might be associated with a temporary lack of an LCV antigen-specific cytotoxic T-cell adaptive immune response. Consistent with the established mechanisms of LCV-related lymphoproliferation, MYC and BCL2L11 gene expression were increased and decreased, respectively. While there was no overt immunosuppression, immunophenotyping revealed the index animal had a relatively low NK cell count, which further decreased by >50% on day 24 of the study. In addition to the temporary lack of adaptive immunity, the low NK cell counts were suggestive of an impaired innate immunity to control the virally-transformed cells and the subsequent unchecked lymphoproliferation. To our knowledge, this is the first report of a Stage V lymphoma with a unique pathogenesis in an otherwise immunocompetent cynomolgus monkey.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins.

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the 'venom-ome' and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 'venom-ome-specific toxins' (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

S100a4 upregulation in Pik3caH1047R;Trp53R270H;MMTV-Cre-driven mammary tumors promotes metastasis.

BACKGROUND: PIK3CA mutations are frequent in human breast cancer. Pik3caH1047R mutant expression in mouse mammary gland promotes tumorigenesis. TP53 mutations co-occur with PIK3CA mutations in human breast cancers. We previously generated a conditionally activatable Pik3caH1047R;MMTV-Cre mouse model and found a few malignant sarcomatoid (spindle cell) carcinomas that had acquired spontaneous dominant-negative Trp53 mutations. METHODS: A Pik3caH1047R;Trp53R270H;MMTV-Cre double mutant mouse breast cancer model was generated. Tumors were characterized by histology, marker analysis, transcriptional profiling, single-cell RNA-seq, and bioinformatics. Cell lines were developed from mutant tumors and used to identify and confirm genes involved in metastasis. RESULTS: We found Pik3caH1047R and Trp53R270H cooperate in driving oncogenesis in mammary glands leading to a shorter latency than either alone. Double mutant mice develop multiple histologically distinct mammary tumors, including adenocarcinoma and sarcomatoid (spindle cell) carcinoma. We found some tumors to be invasive and a few metastasized to the lung and/or the lymph node. Single-cell RNA-seq analysis of the tumors identified epithelial, stromal, myeloid, and T cell groups. Expression analysis of the metastatic tumors identified S100a4 as a top candidate gene associated with metastasis. Metastatic tumors contained a much higher percentage of epithelial-mesenchymal transition (EMT)-signature positive and S100a4-expressing cells. CRISPR/CAS9-mediated knockout of S100a4 in a metastatic tumor-derived cell line disrupted its metastatic potential indicating a role for S100a4 in metastasis. CONCLUSIONS: Pik3caH1047R;Trp53R270H;MMTV-Cre mouse provides a preclinical model to mimic a subtype of human breast cancers that carry both PIK3CA and TP53 mutations. It also allows for understanding the cooperation between the two mutant genes in tumorigenesis. Our model also provides a system to study metastasis and develop therapeutic strategies for PIK3CA/TP53 double-positive cancers. S100a4 found involved in metastasis in this model can be a potential diagnostic and therapeutic target.

Evaluating the quality of the 1000 genomes project data.

BACKGROUND: Data from the 1000 Genomes project is quite often used as a reference for human genomic analysis. However, its accuracy needs to be assessed to understand the quality of predictions made using this reference. We present here an assessment of the genotyping, phasing, and imputation accuracy data in the 1000 Genomes project. We compare the phased haplotype calls from the 1000 Genomes project to experimentally phased haplotypes for 28 of the same individuals sequenced using the 10X Genomics platform. RESULTS: We observe that phasing and imputation for rare variants are unreliable, which likely reflects the limited sample size of the 1000 Genomes project data. Further, it appears that using a population specific reference panel does not improve the accuracy of imputation over using the entire 1000 Genomes data set as a reference panel. We also note that the error rates and trends depend on the choice of definition of error, and hence any error reporting needs to take these definitions into account. CONCLUSIONS: The quality of the 1000 Genomes data needs to be considered while using this database for further studies. This work presents an analysis that can be used for these assessments.

Adaptive adipose tissue stromal plasticity in response to cold stress and antibody-based metabolic therapy.

In response to environmental and nutrient stress, adipose tissues must establish a new homeostatic state. Here we show that cold exposure of obese mice triggers an adaptive tissue remodeling in visceral adipose tissue (VAT) that involves extracellular matrix deposition, angiogenesis, sympathetic innervation, and adipose tissue browning. Obese VAT is predominated by pro-inflammatory M1 macrophages; cold exposure induces an M1-to-M2 shift in macrophage composition and dramatic changes in macrophage gene expression in both M1 and M2 macrophages. Antibody-mediated CSF1R blocking prevented the cold-induced recruitment of adipose tissue M2 macrophages, suggesting the role of CSF1R signaling in the process. These cold-induced effects in obese VAT are phenocopied by an administration of the FGF21-mimetic antibody, consistent with its action to stimulate sympathetic nerves. Collectively, these studies illuminate adaptive visceral adipose tissue plasticity in obese mice in response to cold stress and antibody-based metabolic therapy.

Actionable Activating Oncogenic ERBB2/HER2 Transmembrane and Juxtamembrane Domain Mutations.

Deregulated HER2 is a target of many approved cancer drugs. We analyzed 111,176 patient tumors and identified recurrent mutations in HER2 transmembrane domain (TMD) and juxtamembrane domain (JMD) that include G660D, R678Q, E693K, and Q709L. Using a saturation mutagenesis screen and testing of patient-derived mutations we found several activating TMD and JMD mutations. Structural modeling and analysis showed that the TMD/JMD mutations function by improving the active dimer interface or stabilizing an activating conformation. Further, we found that HER2 G660D employed asymmetric kinase dimerization for activation and signaling. Importantly, anti-HER2 antibodies and small-molecule kinase inhibitors blocked the activity of TMD/JMD mutants. Consistent with this, a G660D germline mutant lung cancer patient showed remarkable clinical response to HER2 blockade.

An Empirical Approach Leveraging Tumorgrafts to Dissect the Tumor Microenvironment in Renal Cell Carcinoma Identifies Missing Link to Prognostic Inflammatory Factors.

By leveraging tumorgraft (patient-derived xenograft) RNA-sequencing data, we developed an empirical approach, DisHet, to dissect the tumor microenvironment (eTME). We found that 65% of previously defined immune signature genes are not abundantly expressed in renal cell carcinoma (RCC) and identified 610 novel immune/stromal transcripts. Using eTME, genomics, pathology, and medical record data involving >1,000 patients, we established an inflamed pan-RCC subtype (IS) enriched for regulatory T cells, natural killer cells, TH1 cells, neutrophils, macrophages, B cells, and CD8+ T cells. IS is enriched for aggressive RCCs, including BAP1-deficient clear-cell and type 2 papillary tumors. The IS subtype correlated with systemic manifestations of inflammation such as thrombocytosis and anemia, which are enigmatic predictors of poor prognosis. Furthermore, IS was a strong predictor of poor survival. Our analyses suggest that tumor cells drive the stromal immune response. These data provide a missing link between tumor cells, the TME, and systemic factors.Significance: We undertook a novel empirical approach to dissect the renal cell carcinoma TME by leveraging tumorgrafts. The dissection and downstream analyses uncovered missing links between tumor cells, the TME, systemic manifestations of inflammation, and poor prognosis. Cancer Discov; 8(9); 1142-55. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 1047.

CRISPR off-target analysis in genetically engineered rats and mice.

Despite widespread use of CRISPR, comprehensive data on the frequency and impact of Cas9-mediated off-targets in modified rodents are limited. Here we present deep-sequencing data from 81 genome-editing projects on mouse and rat genomes at 1,423 predicted off-target sites, 32 of which were confirmed, and show that high-fidelity Cas9 versions reduced off-target mutation rates in vivo. Using whole-genome sequencing data from ten mouse embryos, treated with a single guide RNA (sgRNA), and from their genetic parents, we found 43 off-targets, 30 of which were predicted by an adapted version of GUIDE-seq.

ERK Mutations and Amplification Confer Resistance to ERK-Inhibitor Therapy.

Purpose: MAPK pathway inhibitors targeting BRAF and MEK have shown clinical efficacy in patients with RAF- and/or RAS-mutated tumors. However, acquired resistance to these agents has been an impediment to improved long-term survival in the clinic. In such cases, targeting ERK downstream of BRAF/MEK has been proposed as a potential strategy for overcoming acquired resistance. Preclinical studies suggest that ERK inhibitors are effective at inhibiting BRAF/RAS-mutated tumor growth and overcome BRAF or/and MEK inhibitor resistance. However, as observed with other MAPK pathway inhibitors, treatment with ERK inhibitors is likely to cause resistance in the clinic. Here, we aimed to model the mechanism of resistance to ERK inhibitors.Experimental Design: We tested five structurally different ATP-competitive ERK inhibitors representing three different scaffolds on BRAF/RAS-mutant cancer cell lines of different tissue types to generate resistant lines. We have used in vitro modeling, structural biology, and genomic analysis to understand the development of resistance to ERK inhibitors and the mechanisms leading to it.Results: We have identified mutations in ERK1/2, amplification and overexpression of ERK2, and overexpression of EGFR/ERBB2 as mechanisms of acquired resistance. Structural analysis of ERK showed that specific compounds that induced on-target ERK mutations were impaired in their ability to bind mutant ERK. We show that in addition to MEK inhibitors, ERBB receptor and PI3K/mTOR pathway inhibitors are effective in overcoming ERK-inhibitor resistance.Conclusions: These findings suggest that combination therapy with MEK or ERBB receptor or PI3K/mTOR and ERK inhibitors may be an effective strategy for managing the emergence of resistance in the clinic. Clin Cancer Res; 24(16); 4044-55. ©2018 AACR.

Homozygous KSR1 deletion attenuates morbidity but does not prevent tumor development in a mouse model of RAS-driven pancreatic cancer.

Given the frequency with which MAP kinase signaling is dysregulated in cancer, much effort has been focused on inhibiting RAS signaling for therapeutic benefit. KSR1, a pseudokinase that interacts with RAF, is a potential target; it was originally cloned in screens for suppressors of constitutively active RAS, and its deletion prevents RAS-mediated transformation of mouse embryonic fibroblasts. In this work, we used a genetically engineered mouse model of pancreatic cancer to assess whether KSR1 deletion would influence tumor development in the setting of oncogenic RAS. We found that Ksr1-/- mice on this background had a modest but significant improvement in all-cause morbidity compared to Ksr1+/+ and Ksr1+/- cohorts. Ksr1-/- mice, however, still developed tumors, and precursor pancreatic intraepithelial neoplastic (PanIN) lesions were detected within a similar timeframe compared to Ksr1+/+ mice. No significant differences in pERK expression or in proliferation were noted. RNA sequencing also did not reveal any unique genetic signature in Ksr1-/- tumors. Further studies will be needed to determine whether and in what settings KSR inhibition may be clinically useful.

Silencing of retrotransposons by SETDB1 inhibits the interferon response in acute myeloid leukemia.

A propensity for rewiring genetic and epigenetic regulatory networks, thus enabling sustained cell proliferation, suppression of apoptosis, and the ability to evade the immune system, is vital to cancer cell propagation. An increased understanding of how this is achieved is critical for identifying or improving therapeutic interventions. In this study, using acute myeloid leukemia (AML) human cell lines and a custom CRISPR/Cas9 screening platform, we identify the H3K9 methyltransferase SETDB1 as a novel, negative regulator of innate immunity. SETDB1 is overexpressed in many cancers, and loss of this gene in AML cells triggers desilencing of retrotransposable elements that leads to the production of double-stranded RNAs (dsRNAs). This is coincident with induction of a type I interferon response and apoptosis through the dsRNA-sensing pathway. Collectively, our findings establish a unique gene regulatory axis that cancer cells can exploit to circumvent the immune system.

Single-cell RNA sequencing identifies distinct mouse medial ganglionic eminence cell types.

Many subtypes of cortical interneurons (CINs) are found in adult mouse cortices, but the mechanism generating their diversity remains elusive. We performed single-cell RNA sequencing on the mouse embryonic medial ganglionic eminence (MGE), the major birthplace for CINs, and on MGE-like cells differentiated from embryonic stem cells. Two distinct cell types were identified as proliferating neural progenitors and immature neurons, both of which comprised sub-populations. Although lineage development of MGE progenitors was reconstructed and immature neurons were characterized as GABAergic, cells that might correspond to precursors of different CINs were not identified. A few non-neuronal cell types were detected, including microglia. In vitro MGE-like cells resembled bona fide MGE cells but expressed lower levels of Foxg1 and Epha4. Together, our data provide detailed understanding of the embryonic MGE developmental program and suggest how CINs are specified.